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PREFACE 

The  work  presented  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  under  Program  Element  65807F. 
The  Air  Force  project  manager  was  Mr.  Elton  Thompson  (AEDC/DOTR).  The  results 
presented  were  obtained  by  ARO,  Inc.,  AEDC  Division  (a  Sverdrup  Corporation 
Company),  operating  contractor  for  the  AEDC,  AFSC,  Arnold  Air  Force  Station, 
Tennessee,  under  ARO  Project  Numbers  P32A-C2A  and  P33A-K0A.  Analysis  of  the 
results  was  completed  in  June  1977,  and  the  manuscript  (ARO  Control  No. 
ARO-PWT-TR-77-64)  was  submitted  for  publication  on  September  30,  1977. 
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1.0  INTRODUCTION 

Unsteady  transonic  flow  over  bodies  has  been  treated  with  linear  differential  equations 
for  high  reduced  frequencies  with  major  interest  in  flutter  analysis  (Ref.  1).  Since  the 
frequency  range  of  interest  is  the  low  frequency  range,  where  the  linearized  equation  breaks 
down,  the  nonlinear  equation  must  be  used.  For  slender  bodies  of  revolution  with  a sharp 
nose,  several  investigators  have  developed  approximate  methods  to  solve  the  transonic  small 
disturbance  equation  (Refs.  2 through  4).  For  blunt-nose  bodies  of  revolution,  these 
aforementioned  methods  are  not  applicable.  In  engineering  applications,  a method  that 
can  be  applied  to  both  blunt  and  pointed  bodies  of  revolution  is  most  desirable,  and 
such  a method  is  developed  in  this  report. 

The  method  to  be  described  is  based  on  the  numerical  method  developed  by  South 
and  Jameson  (Ref.  5),  who  solved  the  steady  full  potential  equation  for  ax isym metrical 
flow  over  blunt  and  pointed  bodies  of  revolution  at  transonic  speeds.  Their  method  shows 
good  agreement  with  other  theoretical  methods  and  experiments  at  steady  state  (Ref.  6). 
Thus,  an  extension  of  South  and  Jameson's  scheme  to  calculate  unsteady  transonic  flow 
has  been  pursued. 

The  unsteady  motion  under  consideration  is  assumed  to  be  a small  amplitude  harmonic 
oscillation;  therefore,  the  flow  may  be  divided  into  steady  and  unsteady  components,  with 
the  latter  as  a perturbation  of  the  former.  The  fundamental  equations  for  the  steady  and 
unsteady  parts  were  derived  by  using  a body-normal  coordinate  system  and  are  shown 
to  be  both  elliptic  and  hyperbolic  equations  depending  on  the  local  steady  flow  Mach 
number.  The  numerical  method  of  a rotated  difference  scheme  developed  in  Ref.  5 was 
used  in  solving  both  the  steady  (computer  program  RAXBOD)  and  unsteady  flow  fields. 
Pulsatile  and  pitching  oscillations  were  considered.  Examples  of  the  calculated  results  are 
given  for  both  blunt  and  pointed  bodies  of  revolution  for  quasi-steady  and  unsteady  flow 
conditions. 


2.0  ANALYSIS 


2.1  BASIC  EQUATIONS 

The  success  of  using  a body-normal  coordinate  system  with  rotated  difference  scheme 
in  predicting  surface  pressures  of  blunt  and  pointed  bodies  of  revolution  in  steady  transonic 
flow  (Refs.  5 and  6)  prompts  one  to  use  the  same  approach  for  the  calculation  of  unsteady 
surface  pressure.  The  unsteady  potential  equation  for  an  irrotational  compressible  fluid 
(Ref.  7)  is 
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4 . + 4 t . - v2*  - 4 i .[<?  ■ v)  ?]  <i) 

a2  <?t2  a2  Ot  a2 

where  <p  is  the  velocity  potential,  q is  the  velocity  vector,  a is  the  speed  of  sound,  t 
is  the  time,  and  V and  V2  are  respectively  the  divergent  and  Laplace  operators.  If  Eq. 
(1)  is  applied  to  a body-normal  coordinate  system  (Fig.  1)  for  a body  of  revolution  at 
zero  incidence,  one  obtains 
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For  flows  about  a body  of  revolution  undergoing  the  small  amplitude  oscillation  considered 
here,  a cos  ^ variation  of  the  flow  quantities  in  the  direction  is  assumed.  It  is  also 
assumed  that  the  solution  can  be  expanded  into  a steady  and  an  unsteady  component 
(Ref.  4);  i.e., 
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H ■ I - — n 
ds 

Z ■ J cos  0 ds  - n sin  0 
r - / sin  0 ds  + n cos  0 


Figure  1.  Body-normal  coordinate  system. 


<£  = <£0  (**n)  + [ji  j (s,  n)  cos  V'  eiklJ  R p (4a) 

U = u0  + a*!  costfeikl)R  p (4b) 

V = Vo  + (Vj  cos^  eikt)R  p>  (4c) 

P = P0  + (px  cos  c/r  elkt)R  p (4d) 


where  k is  the  reduced  frequency  of  the  oscillation,  R.P.  signifies  the  real  part  of  a complex 
quantity,  and  the  subscripts  0 and  1 represent  the  quantity  for  steady  and  first-order 
perturbation  flow,  respectively.  By  substituting  Eqs.  (4a  through  d)  into  Eq.  (2)  and 
collecting  terms  of  the  same  order,  one  obtains  the  zeroth-order  equation  for  4>0  as  follows: 


7 


AEDC-TR-77-100 


i 

The  first-order  equation  for  <j>i  is 


where 


A = 2^o\  T sing  _ A (±  - tfP+o  \ 1° 

a2H  r a2H  <?s  \H  5s  / \5sdn/  a2 

0 0 o 


The  quantity  C equals  0 or  1 , depending  on  the  type  of  flow  under  consideration. 
For  axisymmetrical  oscillatory  flow  such  as  the  pulsatile  oscillation,  <j>i  is  independent 
of  0;  therefore,  cos  0 should  be  replaced  by  unity  in  Eq.  (4),  and  the  term  (1/r2  0i) 
should  be  eliminated  from  Eq.  (7);  that  is,  C = 0.  For  a pitching  oscillation,  C = 1. 
Since  Eq.  (7)  is  a complex  equation,  0i  is  a complex  function  and 

01  = 0ir(s,ji)  -*-  i0n  (s,  n)  (10) 

Consequently,  Ui,  Vi,  and  Pi  are  all  complex  functions. 
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2.2  BOUNDARY  CONDITIONS 

The  boundary  condition  is  applied  by  requiring  that  the  velocity  vector  be  tangential 
to  the  surface,  or 

Dl=£!  + uaI  + V(2lL  = 0 fill 

Dt  dt  ds  dn 

and 

B = n - n j (s,  t)  (12) 

Thus,  the  zeroth-order  steady  part  boundary  condition  is 

V = _ sin  0 = 0 (13) 

° 

and  the  first-order  unsteady  part  boundary  condition  is 


d"!  _ Uq  dni  ± d(f>  1 
at  h as  da 


(14) 


Also,  at  infinity  (n  ->  -),  the  velocity  potential  is  required  to  vanish;  that  is,  <t>0  0 

and  <j>i  -*■  0. 

2.3  PRESSURE  COEFFICIENTS 


The  pressure  coefficients  may  be  obtained  by  applying  a perturbation  approach  to 
the  unsteady  Bernoulli  equation,  as  follows: 


y 


1°  = l + izi  M2  [(l  - U*  - v*)]  y"1 

(15) 

p 2 “ L J 

30 

pj_  = _K  (d±i  + ^ + v dJh\ 

p..  p—  \ at  h as  ° an  / 

(16) 

u 

< 

+ 

U 

II 

CJ 

(17a) 

p p0  Pi 

AC  = C„  + iAC 

P1  P1R  P1I 

(17b) 

where 
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and 


(17c) 


(17d) 


(17e) 


3.0  COMPUTATIONAL  SCHEME 

The  numerical  method  employed  for  the  solution  of  the  steady  and  unsteady  flow 
fields  is  the  same  as  that  described  in  Ref.  5.  In  fact,  the  computer  program  RAXBOD 
developed  by  South  at  NASA  Langley  Research  Center  was  used  to  solve  the  steady 
components,  Eqs.  (5)  and  (13).  Figure  2 shows  the  transformation  from  the  physical  plane 
to  the  computational  plane.  The  unsteady  components  are  given  by  Eqs.  (7)  and  (14). 
It  should  be  noted  that  Eq.  (7)  retains  the  same  transonic  flow  characteristics  of  a 
mixed-type,  elliptic/hyperbolic  differential  equation,  depending  on  whether  the  local 
steady  flow  is  subsonic  or  supersonic.  Hence,  the  scheme  given  in  Ref.  5 was  applied 
to  solve  Eqs.  (7)  and  (14).  However,  since  Eqs.  (7)  and  (14)  are  linear,  the  finite  difference 
equations  obtained  are  a set  of  linear  algebraic  equations  and  can  be  solved  by  a direct 
matrix  inversion  instead  of  a relaxation  method,  which  was  used  in  the  solution  of  the 
steady  flow  field  (Ref.  5).  ' 


Computational  Mesh 


Figure  2.  Transformation  from  physical  to  computational  plane. 
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The  details  of  the  application  of  the  rotated  difference  scheme  and  the  relaxation 
method  for  solving  the  steady  flow  field  may  be  found  in  Ref.  5 and  will  not  be  repeated 
herein.  The  implementation  of  the  rotated  difference  scheme  to  the  perturbed  flow,  Eq. 
(7),  and  the  boundary  condition,  Eq.  (14),  as  well  as  the  system  of  algebraic  equations, 
is  given  in  this  section. 

3.1  DIFFERENCE  EQUATIONS 

To  begin  with,  a simple  case  of  a uniform  cylindrical  afterbody  will  be  considered. 
Here,  the  sheared  coordinates  as  described  in  Ref.  5 need  not  be  used.  When  one  applies 
the  rotated  difference  scheme  as  described  in  Ref.  5,  Eq.  (7)  is  rewritten  as  follows: 


H) 


IK. 


I-  I n + C,  1 . folR  + C. 
"'nn  " ~ " 
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where 


-C,  1 - C7  ^1A  + C„  = 0 

6 H ds  1 dn  8 11 


(18b) 
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Similar  equations  may  be  written  for  0nss  and  0nNN.  Now,  the  elliptic  or  hyperbolic 
nature  of  Eqs.  (18a)  and  (18b)  is  determined  by  whether  the  local  value  of  (1  - q^/a2) 
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is  positive  or  negative,  respectively.  Thus,  for  all  terms  in  Eq.  (19a),  if  (1  - I &l)  > 
0,  the  central  difference  should  be  used,  and  if  (1  - q^/a*)  < 0,  the  upwind  difference 
should  be  used.  Furthermore,  the  upwind  quadrant  is  determined  by  the  sign  of  the  local 
V0  component  as  shown  in  Fig.  3.  The  corresponding  different  forms  for  terms  in  Eqs. 
(18a)  and  (18b)  at  the  grid  point  (i,  j)  are  expressed  in  Appendix  A,  and  the  final  difference 
equations  are  given  as  follows: 


For  0 jj 

Supersonic  Points: 


J‘+  1 
j + 2 

i-2  i-1  i 


IfV  > 0 


i 

\ 

i 

V 

1 _ 

f 

■2  i- 

-1  i 

IfV  < 0 


Subsonic  Points: 


j-1 


j 


j+1 


,U . 

t 

Also  for  0 


i-1  i i + 1 

Figure  3.  The  neighboring  points  used  in  the  difference  scheme. 

For  a subsonic  point 


(01  + <*i)  ^lR  + [(06  + a6  - CP  01R  - C6  ^11  ] 

1 — ■ J — 1 1—1,  j 
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1 i+l,j 

+ (0g  - flg)  =0 


(20) 
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For  a supersonic  point  and  V > 0 

>'  1 **  1 H „ . T "l^lH  , , * ^a6  + y 2 “ C4^1K  “ C6^1ll  . 

1—  2 , J 1 — 1 f J — 1 »—  1 . J 
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i . J i.J-rl 

*■  ^6  ^ ] H „ Ka4  + C4^1H  * C6  ^[l  . 

i,j+2  i+l,  j 

f ®9  H = 0 (21) 

1 i+l,j+i 

For  a supersonic  point  and  V < 0 

y,  0]„  + (a2  - yp^iR.  ...  4 l^a6  + ~ C4)^’lH 

i — 2 , j i — 1 , j — 1 

“ C6  ^1l]  . * a7^1U  . . + ye^lR.  , 

i-l, J 1—1  . J—  1 i ,J— 2 

+ [(«2  ~ Vs  + + a3;o1H  + C,  cf>  J | ] 

i.i-l 

+ [(a5  - y3  Cg  + a j + C^)  6j  H - Cg6j[J 

1-J 

+ [(a-  - C:  + a9)  - C;<6,]  =0  (22) 

O 0 7 1 i,J+l 

Similar  equations  can  be  written  for  0n  by  switching  <j> ir  and  <j>u  and  changing  the  sign 
for  C$  and  C7  in  Eqs.  (20)  through  (22);  thus  these  equations  will  not  be  repeated 
herein. 

Special  treatment  is  required  for  the  perturbed  flow  0i  along  the  axis,  as  was  the 
case  for  the  steady  part  4>0.  For  axisymmetric  unsteady  flow  such  as  a pulsatile  oscillation, 
for  example,  s -*■  0, 


and 


0, 


r -»  H0s  + 0,  6 


2, 

2 


d±o 

ds 


0 


(23a) 


U 


(23b) 
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Equation  (18)  reduces  to 


1 /^Vir  \x  /1 

L - Vo\ 

, n ^iR  * 2kVo  d<t>u 

h2  \ dS2  ) \ 

0 

<) 

dn2 

an  .a  <5n 

O 

1 /*Vi,\  + i 

fi  - 

\ *a*i, 

. R ^11  _ 2kVo  **1R 

H2  V ds2  ) ’ 

0 

) dn2 

<9n  a2  dn 

0 

(24a) 

(24b) 


where 


Equation  (24)  should  be  used  along  the  axis,  or  i = 1.  If  the  perturbed  flow  is  not 
axisymmetric,  as  in  the  case  of  pitching  or  heaving  oscillations,  then  because  of  the 
assumption  that  fa  is  antisymmetrical  with  respect  to  i//  (i.e.,  fa  'v  cos  0),  fa  r and 
0ii  must  be  zero  along  the  axis. 


Again,  the  upwind  difference  should  be  applied  to  the  second  term  in  Eqs.  (24) 
when  (1  -V2/a„)  < 0 and  the  central  difference  should  be  used  when  (1  - V20/a.l)  > 0. 
For  other  terms,  the  central  difference  should  be  used.  Hence,  at  each  grid  point,  two 
algebraic  equations  are  obtained  for  fa  r and  fai,  and  these  two  equations  are  coupled 
through  the  terms  with  coefficient  C$  and  C7.  For  points  involving  boundaries,  such  as 
the  body  surface,  the  axis  of  symmetry,  and  infinity,  the  proper  boundary  condition 
should  be  imposed  (to  be  discussed  later).  Hence,  the  total  number  of  algebraic  equations 
to  be  solved  simultaneously  will  be  2 x (I  - 1)  x (J  - 1),  where  I and  J are  the  maximum 
grid  index  along  the  s and  n directions,  respectively. 


3.2  BOUNDARY  CONDITIONS 


The  surface  boundary  condition,  Eq.  (14),  may  be  rewritten  as  follows: 


fom  = Ha  . 

(Jn  H ds 


(25a) 


dn 


k Uj 


(25b) 


where  ni  is  given  by  the  motion  of  the  body  and  will  be  described  later.  In  difference 
form,  Eqs.  (25a)  and  (25b)  become 
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ii. -“£/!!•  - <*,„ 

'■Jt1  S,  \ U 0s  / i.J-l 


(26a) 


4> 


n . 


= - (kn ,)  + <£ 


1'  r S°1I 


i.J+l  S,  1 "i.J-l 

Equations  (26a)  and  (26b)  will  be  used  to  replace  all  terms  at  J+l. 


(26b) 


The  expressions  for  nj  for  three  oscillatory  motions  which  have  engineering 
applications  are  now  described. 


(1)  Pulsatile  oscillation  with  the  amplitude  proportional  to  the  local  radius: 


n , = Sa  (r  cos  0) 


d”  l = S (K  r sin  0 + II  cos  0 sin  0) 

0S  0 


(2)  Pitching  oscillation  with  amplitude  t0  : 


n i = rQ  (<i,  - Z)  cos  0 


= r (d,  - Z)  K sin  0 - cos20 
0s  0 1 


(27a) 

(27b) 

(28a) 

(28b) 


where  dt  is  the  pitch  center. 

(3)  Heaving  oscillation  with  amplitude  h0 : 

n j = Kq  cos  0 


^1  = h K sin  0 
0s  ° 


(29a) 

(29b) 


The  boundary  condition  at  infinity  is  #ir  = 0n  = 0 or  all  values  of  0ir  and  <pu 
at  j = 1 are  zero. 


The  boundary  condition  at  the  axis  for  axisymmetric  flow  becomes 

= ^IR.  , . 
i+l.  i i-l  .j 


(30a) 


^ii  J.1  . = ^ii.  . 

i+i. j i-i  >j 


(30b) 
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For  nonaxisymmetric  flow,  one  obtains 

<$11,  = ©in  and  61r  = 0 (31a) 

i+L..  i-l.J  i.J 

= -^n  , and  tfii.  . = 0 (31b) 

i+l.J  :-l.)  i,j  v ' 

because  of  the  cos  1 \>  variation  of  the  flow  in  the  crossflow  plane.  The  set  of  linear  algebraic 
equations  was  solved  by  the  method  of  factorization. 

4.0  RESULTS  AND  DISCUSSION 

4.1  GENERAL  REMARKS 

A computer  program  to  solve  for  the  perturbed  flow  field  was  coded  for  an  IBM 
370-165  computer.  The  computational  procedure  was  first  to  compute  the  steady  flow 
field  using  the  computer  program  RAXBOD.  The  values  of  0O  in  each  node  and  other 
required  information  were  stored;  then  all  the  coefficients  for  the  algebraic  equation  and 
the  perturbed  boundary  conditions  for  a given  oscillatory  motion  were  computed.  Finally, 
a simultaneous  solution  of  the  algebraic  equations  gave  the  values  of  0i  r and  0i  ] at 
every  nodal  point.  All  perturbed  flow  quantities  could  now  be  obtained  by  differencing 
the  0 o,  0i  r,  and  0u  values. 

Regarding  the  selection  of  the  nodal  distribution  to  give  good  resolution,  two  sets 
of  97  x 33  and  97  x 97  were  frequently  used  for  the  steady  part  calculation.  For  the 
perturbed  unsteady  part,  four  nodal  distributions  were  tried  (i.e. , 25  x 9,  49  x 17,  97 
x 17,  and  97  x 33).  The  results  obtained  from  25  x 9 were  not  satisfactory  in  many 
cases,  particularly  at  a supersonic  free-stream  velocity.  The  results  obtained  from  the  49 
x 17,  97  x 17,  and  97  x 33  matrices  have  good  agreement  in  many  cases.  The  required 
computer  time  is  about  1-1/2  min.  for  49  x 17,  4 min.  for  97  x 17,  and  30  min.  for 
97  x 33.  It  was  decided  that  the  nodal  distribution  of  97  x 17  was  sufficiently  good 
from  the  standpoints  both  of  accuracy  and  economy.  Therefore,  most  results  presented 
in  this  report  are  obtained  using  97  x 17  nodes. 

4.2  PULSATILE  OSCILLATION 

A body  of  revolution  undergoing  oscillatory,  pulsatile  motion  of  the  surface  is  a 
fundamental  axisymmetric,  unsteady  motion.  For  the  perturbation  method  developed  in 
this  report,  the  pulsatile  oscillation  can  also  provide  a direct  comparison  between  the 
quasi-steady  solution  of  the  linear  perturbed  equation  (Eq.  (7))  and  the  steady-state  solution 
of  the  nonlinear  equation  (identical  to  Eq.  (5))  for  the  same  body  geometry.  The  case 
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of  a hemisphere-cylinder  was  chosen  for  comparison.  The  magnitude  of  the  pulsatile 
oscillation  is  assumed  to  be  proportional  to  0.1  r,  where  r is  the  local  radius  of  the  body. 
As  shown  in  Fig.  4,  the  expanded  (Rt  = 1.1)  and  the  shrunken  (Rt  = 0.9)  bodies  are 
an  ellipse-cylinder  with  nose  axis  ratios  of  1.1  and  0.9,  respectively.  The  steady-state 
solution  for  the  surface  pressure  shown  in  solid  curves  in  Figs.  4 (M^  = 0.6)  and  5 (M, 
= 1.2)  for  the  ellipse-cylinders  was  obtained  directly  from  the  computer  program  RAXBOD. 
The  quasi-steady  solution  for  the  same  ellipse-cylinders  can  be  obtained  by  adding  (for 
Rt  = 1.1)  or  subtracting  (for  Rt  = 0.9)  the  perturbed  pressure  from  the  mean  steady-state 
surface  pressure  for  a hemisphere-cylinder  (Rt  = 1.0).  The  quasi-steady  solutions,  as  shown 
by  broken  curves  in  Figs.  4 and  5,  are  seen  to  agree  well  with  the  steady-state  solution. 


Figure  4.  Quasi-steady  surface  pressure  perturbation  for 
hemisphere-cylinder  at  = 0.6. 

4.3  PITCHING  OSCILLATION 

Calculations  for  a pitching  oscillation  were  performed  for  five  body  configurations: 
(a)  4-caliber  ogive-nose-cylinder  with  2/R  = 20;  (b)  2-caliber  ogive-nose-cylinder  with  2/R 
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Figure  5.  Quasi-steady  surface  pressure  perturbation  for 
hemisphere-cylinder  at  M_  = 1 .2. 

= 14;  (c)  2:1  ellipse-nose-cylinder  with  2/R  = 20;  (d)  hemisphere-cylinder  with  C/R  = 
14;  and  (e)  parabolic  arc  nose  with  fi/R  = 5.  The  quasi-steady  (k  = 0)  solutions  for  the 
surface  pressure  will  be  compared  to  steady  experimental  data  for  configurations  a,  c, 
and  d at  zero  and  nonzero  incidence.  The  quasi-steady  aerodynamic  coefficients  and  the 
nonsteady  damping-in-pitch  coefficients  will  be  compared  to  experimental  data  and 
available  theory  for  configurations  b,  d,  and  e. 

4.3.1  Quasi-Steady  Perturbation  Pressure 

4.3.1. 1 The  4-Caliber  Ogive-Nose-Cylinder 

In  Fig.  6,  the  steady-state  solution  for  surface  pressure  (M,  = 0.8  to  1.1)  for  the 
4-caliber-ogive-nose-cylinder  compares  satisfactorily  with  the  experimental  data  (Ref.  8)  at 
zero  incidence.  Thus,  the  steady-part  solution,  which  was  used  as  the  input  for  the  unsteady 
part  calculation,  is  valid.  The  quasi-steady  solution  for  the  perturbed  surface 
pressure  is  shown  in  Fig.  7 for  the  leeward  plane  of  symmetry  (all  discussion  in  this 
report  is  referred  to  the  leeward  plane  of  symmetry).  Since  the  perturbed  flow  field  is 
linear  with  respect  to  the  incidence,  the  magnitude  of  A Cp/a  is  plotted.  In  Ref.  8,  two 
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C 


P 


Figure  6.  Comparison  of  steady-state  surface  pressure  between  theory 
and  experiments  for  a 4-caliber  ogive-nose-cylinder. 


4-Caliber  Ogive  Nose 
Theory—  Exp  ° J S 


■I.01 ■ ■ ■ * 1 1 ■ • 1 ■ 1 
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Z/R  Z/R 


Figure  7.  Comparison  of  quasi-steady  perturbed  surface 
pressure  between  theory  and  experiments  for  a 
4-caliber  ogive-nose-cylinder. 
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sets  of  experimental  data  were  obtained  for  a = ±4  and  ±8  deg.  To  obtain  the  average 
perturbed  pressure  slope  for  comparison,  the  following  formula  was  used.: 


AC  (C  „ - C ) 
p __  v pa  p —a  ' 

a 2a 


(32) 


Figure  7 shows  that  the  agreement  between  theory  and  experiment  is  excellent.  That  the 
data  for  a = ±8  and  ±4  deg  agree  well,  except  at  = 1.1,  indicates  the  linearity  of 
the  perturbed  flow  up  to  8 deg 


4.3.1.2  The  2:1  Ellipse-Nose-Cylinder 

In  Fig  8,  the  steady-state  solutions  for  surface  pressure  (Mm  = 0.6  to  1.2)  for  the 
2:1  ellipse-nose-cylinder  again  compare  satisfactorily  with  the  experimental  data  (Ref.  8) 


Figure  8.  Comparison  of  steady-state  surface  pressure 
between  theory  and  experiments  for  a 2:1 
ellipse-nose-cylinder. 
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at  zero  incidence.  At  incidence,  the  experimental  data  were  obtained  for  = 0.9  to 
1.2  and  a = ±8  deg  only.  The  comparison  between  the  quasi-steady  theory  and  experiments 
for  the  average  perturbed  surface  pressure  is  given  in  Fig.  9a.  The  agreement  is  generally 
good.  However,  the  theoretical  data  show  a spike  and  dip  for  Mb  = 0.9  and  a hump 
for  > 1.0  aft  of  the  juncture  of  the  ellipse  and  the  cylinder.  These  phenomena  are 
not  shown  by  the  experimental  data.  The  cause  may  be  attributed  to  the  discontinuity 
in  the  curvature  of  the  body  geometry  in  the  inviscid  calculation  of  the  perturbed  flow 
field.  In  order  to  show  the  effect  of  curvature  discontinuity,  plots  of  the  experimental 


Theory  — 
Exp  °ta 


2: 1EI I ipse  Nose 
□ g aegfCZZJ 


a.  Average  perturbed  surface  pressure 
Figure  9.  Comparison  of  quasi-steady  perturbed  surface 

pressure  between  theory  and  experiments  for  a 2:1 
ell  i pse-nose-cyl  inder . 

leeward  side  perturbation  pressure  (CPa=8°  - CPa  = Q)  and  windward  side  perturbation 
pressure  (CPa=0  - CPa=8»)  are  presented  in  Fig.  9b.  Interestingly  enough,  the  experimental 
data  for  the  windward  side  do  show  the  expected  type  of  peak  and  dip  for  Mm  = 0.9 
and  the  expected  hump  for  M.  > 1.0  approximately  in  the  locations  given  by  the  theory, 
but  to  a far  lesser  degree  than  indicated.  The  experimental  data  in  the  leeward  side  behave 
in  a significantly  different  manner  from  those  of  the  windward  side  in  2 < z/R  < 2.5 
for  from  0.9  to  1.1,  indicating  the  effect  of  viscosity  and  hence  the  deviation  from 
the  assumption  of  a cos  i//  flow  field  at  a = 8 deg. 
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o o • -8  deg  (Windward  Side) 
Theory  Exp  a ~ 8 deg  (leeward  Side) 


b.  Windward  and  leeward  side  perturbed  surface  pressure 
Figure  9.  Concluded. 


4.3.1. 3 The  Hemisphere-Cylinder 

The  third  configuration  used  for  surface  pressure  comparison  is  a hemisphere-cylinder. 
A good  agreement  in  the  surface  pressure  between  theory  and  experiment  at  zero  incidence 
when  there  is  no  flow  separation  has  been  shown  in  Ref.  9,  but  flow  separation  does 
occur  at  = 0.7  to  0.9  at  zero  incidence.  To  evaluate  the  quasi-steady  theory  at  incidence, 
the  Ref.  9 experimental  data  for  hemisphere-cylinders  at  a = ±5  deg  were  used.  (In  Ref. 

9,  <j>  = 0 and  180  deg,  which  corresponds  to  a = 5 and  -5  deg,  used  herein.)  In  Fig. 

10,  the  leeward  and  windward  side  perturbation  pressures  are  compared  between  theory 
and  experiment  for  = 0.6,  0.7,  0.8,  0.9,  1.0,  and  1.2.  In  Fig.  11,  the  shadowgraphs 
of  the  flow  field  at  a = 5 deg  are  shown.  At  = 0.6,  Fig.  1 1 shows  that  the  flow 
is  attached  except  for  a possible  minor  separation  near  the  juncture  of  the  hemisphere 
and  the  cylinder.  This  separation  is  also  indicated  by  the  experimental  data  of  the 
perturbation  pressure  on  the  leeward  and  windward  sides.  The  theory  seems  to  slightly 
overpredict  the  perturbation  pressure.  For  = 0.7  to  0.9,  Fig.  1 1 shows  that  the  leeward 
flow  separates.  On  the  windward  side,  a small  shock  system  is  indicated  for  = 0.7; 
for  Mb  = 0.8,  the  shock  system  resembles  the  case  of  a = 0,  which  indicates  the  existence 
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of  a separation  bubble  (see  Ref.  9).  For  = 0.9,  an  embedded  normal  shock  is  shown. 
For  the  described  flow  field,  one  would  not  expect  good  agreement  between  the  theory 
and  the  experiments  because  of  viscous  effects.  As  shown  in  Fig.  10,  the  agreement  is 


Figure  10.  Comparison  of  quasi-steady  perturbed  surface 
pressure  between  theory  and  experiments  for  a 
hemisphere-cy  I inder. 
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poor  for  these  cases  just  discussed.  It  is  noted  that  there  are  spikes  in  the  theoretical 
curves  for  = 0.7  to  0.9  at  the  region  near  the  zero  incidence  shock  location.  The 
calculated  values  at  the  spikes  are  unrealistic  because  of  the  large  gradients  in  the  flow 
quantities  near  the  shock  (see  Section  4.3.4  for  modification).  For  Mm  = 1.0  and  1.2, 
Fig.  1 1 shows  that  the  flow  is  again  well  attached,  which  is  also  indicated  by  the 
experimental  data  shown  in  Fig.  10.  Thus,  for  Mm=  1.0  and  above,  the  comparison  between 
theory  and  experiments  is  good  except  for  the  hump  resulting  from  the  curvature 
discontinuity  aft  of  the  junction  of  the  hemisphere  and  the  cylinder,  where  the  theory 
shows  a larger  value  again. 

o a ■ -5  deg  (Windward  Side) 

Theory  txp  g a „ 5 deg  (Leeward  Side) 


0.2 


Figure  10.  Concluded. 
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Mqq  = 0. 6 Mm  * 0.7  Mqj-0.8 


© 

Flow  Separation 

® 

Embedded  Shock  System 

© 

Embedded  Normal  Shock 

© 

Recompression  Shock 

© 

Intersection  of  Shock  and  Tunnel  Side  Walls 

Bow  Shock 

© 

Disturbances  Reflected  from  Wall  Due  to  Bow  Shock 

Figure  11. 

Shadowgraphs  of  hemisphere-cylinder  at  5 deg 

4.3.2  Shock  Excursion  at  Incidence 

The  surface  pressure  obtained  by  the  perturbation  theory  and  the  experimental  data 
are  presented  in  Fig.  12  for  the  ellipse-cylinder  at  = 0.9  and  a = 0 and  ±8  deg  (no 
serious  flow  separation  occurs  on  the  leeward  side  of  this  case).  At  a = 0 the  shock 
occurs  at  approximately  Z/R  = 2.3,  where  the  pressure  increases  sharply.  It  is  noted  that 
the  experimental  data  for  a = ±8  deg  cross  over  the  data  for  a = 0.  This  is  an  indication 
of  the  moving  of  the  shock  as  the  body  pitches  to  an  angle.  For  example,  if  the  body 
is  at  positive  incidence  (leeward  side),  the  shock  will  move  forward  from  0 to  1 as  shown 
in  the  sketch  in  Fig.  1 2.  The  surface  pressure  in  between  the  points  0 and  1 will  increase 
due  to  the  forward  moving  of  the  shock,  whereas  upstream  of  1 the  pressure  will  decrease 
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Figure  12.  Comparison  of  surface  pressure  between  theory  and 
experiments  for  a 2:1  ellipse-nose-cylinder  at 
incidence  and  M,,,  = 0.9. 

due  to  the  incidence.  For  negative  incidence  (windward  side)  the  shock  moves  to  2;  then 
the  pressure  upstream  of  0 increases,  but  the  pressure  in  between  0 and  2 should  decrease 
due  to  the  aft  movement  of  the  shock  from  0 to  2.  Thus,  when  the  data  at  incidence 
cross  over  at  zero  incidence,  this  point  of  crossover  may  be  interpreted  as  the  location 
of  the  shock  at  incidence.  The  excursion  of  the  shock  as  read  by  this  method  is  shown 
in  Fig.  12.  Since  the  distance  between  points  1 and  2 is  about  0.2R  and  covers  a total 
incidence  of  16  deg,  the  shock  excursion  is  relatively  small  for  bodies  of  revolution  as 
compared  to  a 2-D  airfoil  for  which  a 10-percent  chord  length  excursion  over  a few  degrees 
of  incidence  is  common.  In  addition,  the  theoretical  result  shows  a point  of  crossover 
in  the  surface  pressure  near  the  middle  of  points  1 and  2 and  has  the  same  trend  as 
shown  by  the  experimental  data.  Therefore,  the  presumed  fixed  shock  location  in  the 
development  of  the  perturbation  theory  for  bodies  or  revolution  is  not  a serious  deficiency 
when  the  shock  is  not  too  strong.  If,  however,  the  normal  shock  is  strong  (for  increasing 
bluntness  of  the  nose),  as  for  the  case  of  a hemisphere-cylinder,  the  flow  in  the  leeward 
side  is  generally  separated,  and  the  inviscid  calculation  is  invalidated.  For  the  windward 
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side,  however,  the  shock  excursion  is  not  large  (about  0.2R),  as  shown  in  Fig.  10  for 
= 0.9. 

4.3.3  Quasi-Steady  Aerodynamic  Coefficients 

The  integrated  aerodynamic  coefficients  for  the  2-caliber  ogive-nose-cylinder  and  the 
hemisphere-cylinder  are  shown  in  Figs.  13  and  14,  respectively.  The  experimental  force 
data  were  obtained  from  Refs.  9 and  10.  The  theoretical  value  of  Ca,f  is  the  same  as 
at  a = 0 and  is  generally  lower  than  the  experiments.  However,  experimental  force  data 

Fores  Data,  Ref.  10  ^ ^ X/R  ■ 20 

Inviscid,  X/R  ■ 14 

Theory  Viscous  Effect  Added 


Mco 


Figure  13.  Comparison  of  static  stability  derivatives 

between  theory  and  experiments  for  a 2-caliber 
ogive-nose-cylinder. 
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Force  Data,  Ref.  10  | " ?• 

\±IIR » 18 

Pressure  Data,  Ref.  9 1 & ilR  ■ 20 

Theory  Invlscld,  IlR  - 14 

Viscous  Effect  Added 


Figure  14.  Comparison  of  static  stability  derivatives 
between  theory  and  experiments  for  a 
hemisphere-cylinder. 

for  the  shorter  body  and  the  integration  of  surface  pressure  measurements  for  the 
hemisphere-cylinder  (Fig.  14)  do  show  better  agreement  with  the  theory.  This  indicates 
that  the  difference  is  due  to  the  skin  friction,  which  was  not  accounted  for  in  the  theory. 
When  the  viscous  effect  is  added,  using  a skin  friction  corresponding  to  the  experimental 
Reynolds  number  for  the  long  body  only,  the  agreement  between  theory  and  experiments 
for  CA  ,f  is  excellent.  The  CN  a value  is  nearly  independent  of  the  body  length  and  is 
shown  to  be  in  good  agreement  with  the  theory  for  the  ogive-cylinder.  For  the 
hemisphere-cylinder,  the  hump  and  dip  in  the  Cn  a data,  given  by  the  theory  as  between 
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Mach  numbers  from  0.6  to  0.9,  were  not  shown  by  the  experimental  data;  but  for  Mm 
>0.9  the  agreement  is  again  good.  The  aerodynamic  moment  center  (neutral  point) 
measured  from  the  body  midpoint  (positive  forward)  is  compared  between  theory  and 
experiments  in  Figs.  13  and  14;  the  agreements  are  only  fairly  good. 

4.3.4  Modification  at  th6  Region  of  Embedded  Shock 

As  described  in  Section  4.3.1,  unrealistic  spike  was  found  at  the  region  of  embedded 
shock.  South  and  Jameson  used  the  criterion  of  minimum  Laplacian  (V20o)  to  define 
the  location  of  shock  and  smooth  the  flow  field  immediately  upstream  and  downstream 
of  the  shock  by  one-sided  extrapolation.  In  this  section,  the  same  one-sided  extrapolation 
at  the  embedded  shock  is  used  for  the  quasi-steady  flow  solution,  and  the  effect  of  the 
embedded  shock  on  the  solution  of  the  perturbed  flow  field  is  examined.  The  cases 
examined  are  the  2:1  ellipse-nose-cylinder  at  Mm  = 0.9  and  the  hemisphere-cylinder  at 
M,  = 0.7,  0.8,  and  0.9.  Significant  changes  are  found  in  the  pressure  distribution  near 
the  embedded  shock,  as  the  following  paragraph  describes. 

Theoretical  results  with  and  without  modification  are  shown  in  Fig.  1 5 by  solid  and 
dashed  curves,  respectively.  For  the  ellipse-nose-cylinder  (Fig.  15a),  the  spike  value  is 


Z/R 

a.  Ellipse-cylinder,  = 0.9 

Figure  15.  Comparison  of  theoretical  results  with  and  without 
modification  and  experimental  data. 
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Theory  Exp 


ZIR 

b.  Hemisphere-cylinder 
Figure  15.  Concluded. 


reduced  from  -4.31  to  -1.8,  approximately,  and  in  the  downstream  of  the  shock,  the 
perturbed  pressure  does  not  become  positive,  as  shown  by  the  previous  result  and  the 
experimental  data.  Now,  the  theoretical  result  for  M,  = 0.9  looks  just  like  the  other  cases 
for  > 1.0.  For  the  hemisphere-cylinder  at  = 0.7,  Fig.  15b,  the  spike  value  reduces 
from  -0.84  to  -0.4.  The  new  results  of  surface  pressure  change  significantly  after 
modification  at  M„  = 0.8  and  0.9.  At  M„  = 0.8,  the  spike  value  of  ACP  reduces  from 
1.05  to  0.2  and  then  monotomically  approaches  zero  without  the  second  spike  in  the 
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-ACp  side;  and  at  = 0.9,  the  spike  at  the  -ACP  side  disappears  entirely  while  the 
positive  side  reduces  from  0.64  to  0.24.  The  comparison  between  theory  (after 
modification)  and  experiment  is  fair  in  the  windward  side  except  at  M,,,  = 0.8  where  flow 
separation  prevails  even  at  the  windward  side.  From  the  above  comparison,  it  is  clear 
that  the  perturbed  surface  pressure  obtained  with  modification  at  the  embedded  shock 
provides  better  solution.  It  should  be  pointed  out  that,  although  the  present  results  differ 
significantly  from  the  previous  results  in  surface  pressure  distribution,  the  integrated 
aerodynamic  coefficients  change  only  slightly.  Hence,  the  results  reported  in  Section  4.3.3 
remain  unchanged. 

4.3.5  Comparison  with  Nonlinear  Theory 

Inviscid  nonlinear  computations  made  by  solving  the  Euler's  equations  are  available 
for  a hemisphere-cylinder  at  a = 5 deg  and  Mit=  1.2  (Ref.  9)  and  = 0.9  (Ref.  11). 
In  Fig.  16a,  the  perturbed  surface  pressure  distributions  are  compared  between  the 


Present  Non-Linear  Theory,  Ref.  9 and  11 

Windward  Side 

Leeward  Side 


a.  Quasi -steady  perturbation  pressure 
Figure  16.  Comparison  of  perturbation  theory  with  nonlinear 
(Euler's  equations)  theory  for  a hemisphere-cylinder 
at  a = 5 deg  and  M„  = 0.9  and  1.2. 
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nonlinear  theory  and  the  present  method.  Since  in  Ref.  11  only  the  Cp  values  at  \j/  - 
0 deg,  90  deg,  and  180  deg  were  given,  the  Cp  values  at  ^ = 90  deg  were  used  to  replace 
Cp  values  at  a = 0 to  obtain  the  perturbed  surface  pressure.  That  is, 


Windward  side 

ACp 

= cn  - c 

p 

P^=180°  Ptf/=  90 

Leeward  side 

AC 

= C - c 

p 

pi£  = 90°  PiZr=0 

The  ACP  curves  between  present  calculation  and  the  nonlinear  theory  agree  well  in  the 
trend.  It  is  also  noted  that  at  a = 5 deg  the  assumption  of  cos  \jj  variation  of  flow  variables 
in  the  present  theory  is  perhaps  invalid  at  the  nose  for  = 1.2  and  at  1.5  < z/R  < 
3-0  for  = 0.9. 


b.  Quasi-steady  pressure 
Figure  16.  Concluded. 


A plot  of  the  total  pressure  Cp  for  a hemisphere-cylinder  at  a = 5 deg  and  it  = 
0,  90,  and  180  deg  is  shown  in  Fig.  16b  for  = 0.9  and  1.2.  The  agreement  between 
present  theory  and  the  nonlinear  theory  is  good  except  near  the  embedded  shock  region 
for  M„  = 0.9  and  aft  of  the  juncture  of  the  hemisphere  and  the  cylinder  for  M„,  = 1.2. 
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The  experimental  data  at  = 0 and  90  deg  for  = 0.9  clearly  show  that  the  flow 
is  separated  in  these  two  meridian  planes. 

4.3.6  Unsteady  Calculation 

Unsteady  calculations  were  first  performed  for  a parabolic  arc  nose  at  a reduced 
frequency  of  k = 0.006  and  = 1.0.  This  case  was  computed  by  Liu  and  Ruo  (Refs. 
2 and  12)  by  extending  the  parabolic  method  of  Oswatisch  and  Keune  (Ref.  13).  The 
comparison  of  the  stability  derivatives  of  the  total  damping-in-pitch  from  the  present  theory 
and  the  theory  of  Liu  and  Ruo  as  well  as  the  linearized  theory  of  Landahl  (Fig.  12  of 
Ref.  12)  is  shown  in  Fig.  17  for  a variation  in  pitch  center.  The  agreement  between  the 


Figure  17.  Comparison  of  theoretical  stability  derivatives 
for  a parabolic  arc  nose. 
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Figure  17.  Concluded. 

present  theory  and  that  of  Liu  and  Ruo  is  good  except  for  dt  0,  where  the  present 
theory  gives  a lower  value.  It  should  be  noted  that  Landahl's  linearized  theory  is  invalid 
for  such  a low  reduced  frequency.  Liu's  theory  indicates  a reduction  of  the  total 
damping-in-pitch  coefficient  when  the  nonlinearity  of  the  differential  equation  is 
approximately  accounted  for.  The  present  theory  indicates  that  if  the  nonlinearity  is  fully 
accounted  for,  the  total  damping-in-pitch  coefficient  is  even  lower.  No  experimental  data 
are  available  for  comparison.  Other  stability  derivatives  such  as  CN<I  + CNq,  CNa,  and 
Cm  q calculated  by  the  present  method  as  a function  of  dt  are  also  shown  in  Fig.  1 7. 
Again,  no  data  are  available  for  comparison  for  these  quantities. 

No  transonic  unsteady  surface  pressure  measurements  for  a body  of  revolution  have 
been  reported  in  the  literature.  Therefore,  a comparison  of  the  unsteady  surface  pressure 
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is  not  available.  However,  a stability  derivative  for  total  damping-in-pitch,  Cm  ^ + On  , 
has  been  reported  in  Refs.  14  and  15  for  a 2-caliber  ogive-nose-cylinder  and  a 
hemisphere-cylinder.  In  these  measurements,  the  models  were  pitched  at  a constant 
frequency  of  36  Hz;  hence  the  reduced  frequency  changes  with  the  free-stream  Mach 
number  as  follows: 

= 0.6  0.7  0.8  0.9  1.0  1.1  1.2 

k = 0.0165  0.0145  0.0124  0.0112  0.0099  0.0095  0.0088 

Calculations  were  made  using  the  corresponding  reduced  frequency  for  each  Mach  number. 
Figure  18  shows  the  stability  derivatives  for  total  damping-in-pitch  for  the  ogive-cylinder. 
The  Cm.  + Cmq  value  given  by  the  theory  is  generally  higher  than  the  experimental 
data,  particularly  for  the  lower  Mach  number.  No  data  are  available  for  comparison  for 
cNi  + CNq. 


Theory  M iTR 

O Exp  tyR  ■ 7 — A 


Figure  18.  Comparison  of  dynamic  stability  derivatives  between 

theory  and  experiments  for  a 2-caliber  ogive-nose-cylinder. 

For  the  hemisphere-cylinder,  the  comparison  between  theory  and  experiments  is 
shown  in  Fig.  19  for  two  pitch  centers.  The  experimental  data  for  dt/R  = 7 were  obtained 
from  Ref.  13  for  the  full  span  model;  data  for  dt/R  = 5 and  6 are  from  Ref.  14.  The 
agreement  between  theory  and  experiment  is  fair.  The  theoretical  calculations  were  made 
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for  = 0.6,  0.7,  0.8,  0.9,  1.0,  and  1.2  to  determine  the  curves,  which  at  first  look 

do  not  show  the  trend  given  by  the  experimental  data  for  dt/R  = 5 and  7.  However, 
the  experiments  for  dt/R  = 7 gave  a positive  value  of  Cm^  + Cmq  (unstable)  at  M..  = 
0.9,  which  was  not  predicted  by  the  theory,  and  the  set  of  data  for  dt/R  = 5 also  indicate 
a drop  of  Cm  ■ + Cm  at  = 1 .4.  When  these  characteristics  of  the  experimental  data 
are  considered,  the  theoretical  curves  do  not  seem  unreasonable.  In  addition,  the 
experimental  data  for  the  dt/R  = 6 set  behaved  differently  from  the  other  sets.  Hence, 
perhaps  one  needs  more  experimental  data  in  the  Mach  number  range  shown  in  Fig.  19 
in  order  to  fully  evaluate  the  theoretical  results.  As  pointed  out  previously,  for  a 
hemisphere-cylinder  in  the  transonic  flow  regime,  the  flow  field  is  complicated  by  the 
viscous  and  inviscid  interaction.  Therefore,  inviscid  theory  may  not  be  sufficient  to  predict 
the  dynamic  stability  derivatives.  Theoretical  values  for  Cn^  + CNq  are  also  plotted  in 
Fig.  19. 


Theory  Exp  df/R 


Figure  19.  Comparison  of  dynamic  stability  derivatives  between 
theory  and  experiments  for  a hemisphere-cylinder. 
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5.0  CONCLUDING  REMARKS 

A perturbation  theory  and  numerical  solution,  based  on  the  unsteady  full  potential 
equation,  were  developed  for  unsteady  transonic  flow  about  blunt  and  pointed  bodies 
of  revolution  undergoing  harmonic  oscillations.  The  oscillatory  motion  was  considered  to 
be  a small  perturbation  of  the  nonlinear  steady  flow.  The  coupled  equations  for  the  steady 
and  unsteady  potentials  were  solved  numerically  using  a rotating  difference  scheme. 
Calculations  were  performed  for  ogive-cylinder,  ellipse-cylinder,  hemisphere-cylinder,  and 
parabolic  arc  nose  shapes  undergoing  pulsatile  and  pitching  oscillations.  Comparisons  of 
the  present  theory  with  available  experiments  and  other  theories  are  presented  for  the 
surface  pressure  distribution  at  quasi-steady  state  and  for  the  static  and  dynamic  stability 
derivatives.  The  agreement  is  good  for  the  quasi-steady  cases. 
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APPENDIX  A 

DIFFERENCE  FORM  OF  TERMS  IN  EQ.  {18) 

The  difference  form  of  each  term  involved  in  arriving  at  Eqs.  (20),  (21),  and  (22)  is 
given  below. 

I.  For  the  d2/3S2  term 


A.  If  (1  - q2/a2)  > 0 


where 


+ ^lRi>w  + ^3^1Ri+1>h1  + 04* 


i+l.j 


+ 09^1  R. 


i+ 1 , j+ 1 


(A-l) 


01  = “03  “ “07  _ 09  “ a2^ 

02  = a3  «!j_LQ 
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where  the  terms  At  and  Af  are  the  mesh  spacing  in  the  computational  plane,  Fig.  2.  The 
terms  f and  g are  the  stretching  function  in  the  s and  n directions.  The  expression  H is 
defined  in  Fig.  1,  and  the  terms  i and  j are  the  grid  points  along  the  s and  n directions, 
respectively. 


B.  If  (1  - q2/a|)  < 0 
1.  When  V0  > 0 
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where 
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2.  When  VQ  < 0 or  on  the  boundary,  j = J (see  Fig.  3) 

Q . - y[  ^ir.  , . + y 2 ^ir.  . . + y£  0ir.  . 

\ ssy  i,j  1-2,]  1-l.j  i,j 


+ ^ 4 ^lR.  y 5 ^1R.  . * ^6  ^1R.  . , 

i-I, j-1  i , j— 1 i.j-2 

where 


y[  = y l 


II.  For  the  32/3N2  term 

This  term  is  independent  of  the  local  flow  condition,  or 


(A-3) 
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where 
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t 

The  rest  of  the  terms  in  Eq.  (18)  may  be  written  as  follows: 


c4  . ^ir  _ cr 

H 3s  4 

■ ,4ihw,i) 

(A-5) 

C-  = C'  ( 

3 dn  5 V 

(A-6) 

Ce  _ c; 

H ds  6 

" ’ViJ 

(A-7) 

c tylf  _ c . 

7 dn  7 

(A-8) 

C8 

l = C8  ^lR. 

l.J 

(A-9) 
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NOMENCLATURE 


a 

Speed  of  sound 

Ai,Bj,Cj 

Coefficients 

ai.bi.q 

Coefficients 

B 

Equation  for  the  boundary,  Eq.  (12) 

Ca.f 

Forebody  axial-force  coefficient,  FA  /q 

Cn  aiCm  o 

Stability  derivative  caused  by  angle  of  attack 

CNi  + CNq) 

+ Cmq 

Stability  derivative  caused  by  total  damping-in-pitch 

cp 

Pressure  coefficient,  (p  - p„)/q 

ACp 

Perturbed  pressure  coefficient,  Pi/q 

ACp 

Average  perturbed  pressure  coefficient,  Eq.  (37) 

dt 

Pitching  center  from  the  nosetip 

Fa 

Axial  force  on  forebody 

H 

(1  - n ds/dfi) 

h0 

Amplitude  of  harmonically  heaving  oscillation 

I,  J 

Maximum  grid  index  along  the  s and  n directions, 
respectively 

K 

Curvature  of  body  meridian  plane,  -d0/ds  . 

k 

Reduced  frequency,  ojR/IL 

i 

Total  length  of  the  body 

M 

Mach  number 

n,s 

Body  normal  and  tangent  coordinates 
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ni 

P 

Ap 

Q 

q 

q 

R. 


Rt 


r,z 

t 

U,V,W 

X,Y,Z 


Xn.p 


a 


OiA,7i 

7 

80 


$,T 

e 

9 

9 


To 


Perturbed  body  boundary  due  to  motion 
Pressure 

Perturbed  pressure 

(1  - q?/a g) 

Dynamic  pressure,  1/2 
Velocity  vector 

Radius  of  the  cylindrical  body 

Radius  of  the  expanded  or  shrunken  cylindrical  body 

Radial  and  axial  coordinates 
Time 

Axial,  vertical,  and  circumferential  velocity  components,  respectively 
Cartesian  coordinates 

Aerodynamic  moment  center  or  neutral  point  of  pitching  moment,  positive 

forward  ^f  body  mid-point,  (CM  /C„. ) . , 

■ a mid-point 

Angle  of  attack 
Coefficients 
Gas  constant 

Amplitude  of  harmonically  pulsatile  oscillation 
Computational  coordinates 
Slope  of  body  geometry  in  the  meridian  plane 
Dynamic  viscosity  of  the  fluid 
Density  of  the  fluid 

Amplitude  of  harmonically  pitching  oscillation 
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<t>  Velocity  potential 

Circumferential  angle  starting  from  lee-side  plane  of  symmetry 
co  Oscillation  frequency,  cycle/sec 

SUBSCRIPTS 

0. 1  Steady  and  perturbed  unsteady  components,  respectively 

ij  Index  along  the  body  and  normal  coordinates,  respectively 

1, R  Imaginary  and  real  parts,  respectively 

«•  Free  stream 
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